Controlling extended systems with spatially filtered, time-delayed feedback 
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We investigate a control technique for spatially extended systems combining spatial filtering with 
a previously studied form of time-delay feedback. The scheme is naturally suited to real-time control 
of optical systems. We apply the control scheme to a model of a transversely extended semiconductor 
laser in which a desirable, coherent traveling wave state exists, but is a member of a nowhere stable 
family. Our scheme stabilizes this state, and directs the system towards it from realistic, distant and 
noisy initial conditions. As confirmed by numerical simulation, a linear stability analysis about the 
controlled state accurately predicts when the scheme is successful, and illustrates some key features 
of the control including the individual merit of, and interplay between, the spatial and temporal 
degrees of freedom in the control. 
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I. INTRODUCTION 

Nonlinear dynamical systems often possess periodic or- 
bits that have desirable properties but are unstable. The 
problem of applying small perturbations to the system in 
such a way as to produce stable periodic behavior has re- 
ceived much attention recently. |Q] This paper addresses 
the control problem as it arises in a specific context: the 
stabilization of unstable traveling wave states of spatially 
extended systems. Though such states have a particu- 
larly simple structure, the control problem is nontrivial. 

The general method introduced below may be applica- 
ble to a wide variety of physical systems, but an entirely 
general analysis of it is beyond the scope of this work. 
We have chosen to investigate in detail its application 
to two sets of model equations describing the dynam- 
ics of wide aperture semiconductor lasers. Our results 
demonstrate that unstable traveling wave states can be 
effectively controlled in these systems and therefore have 
implications both for the general theory of control of spa- 
tially extended systems and for the design of semiconduc- 
tor lasers. 

For the purposes of this paper, controlling a system 
means providing feedback that locks the system to one 
member of a possibly infinite family of unstable peri- 
odic orbits present in that system, thereby choosing a 
desired state from a large variety of possibilities. The 
technological goal is to produce a desirable behavior in 
a system by applying carefully chosen feedback that di- 
rects the system to the goal state and keeps it there. For 
many applications, it is desirable to design the feedback 
such that the magnitude of the control signal decreases 
as the system approaches the desired state, and, in the 
absence of noise, vanishes when the controlled behavior 
is a dynamical state of the uncontrolled system. It is also 
worthwhile to consider "controlling" a state which is only 
approximately a true orbit of the system. An important 
example is the situation where a stress is applied to a 
large but finite transverse portion of a system. Useful 
results for the dynamics of the active region may still 
be obtained by using feedback designed for the traveling 



wave solution of the infinite, idealized system. In this 
case one might expect the feedback to become small, but 
not completely vanish. 

For a system with an accessible dynamical field A{x, t), 
our control signal eA(x, t) is derived from an infinite sum 
of signals delayed in time by integer multiples of the pe- 
riod of the state that is to be stabilized: 



eyi(x,i) 



(1) 



A(x, t) - (1 - i?) ^ i?"-ii„(x, t 



where 7 is the gain of the feedback, t is the period of 
the target state, and < i? < 1 determines the rela- 
tive weight given to states farther in the past. The field 
Anix,t) ^ T-\p'\k~k^)T{A[x,t)]\ is the spatially fil- 
tered version of A. Here /() is a filtering function ap- 
plied in Fourier space and T\A\ is the spatial Fourier 
transform of A. The precise manner in which the spatial 
filtering is included may vary; we have made a choice that 
corresponds directly to an experimental arrangement de- 
scribed below. We take /(g) to be peaked around zero so 
that contributions to A from wave numbers other than 
the desired wave number, k^, are suppressed in A. We 
also take /(O) = 1 so that the feedback term vanishes 
identically when the system is in a pure state of wave 
number /cc that is an unstable orbit of the uncontrolled 
system. Eq. |l] represents the enhancement of time-delay 
feedback of the form analyzed in Ref. [|| with spatial 
filtering of the type introduced in Ref. H . 

Control based on Eq. |^ is especially well-suited to spa- 
tially extended states with a structure dominated by one 
Fourier mode. Feedback occurs whenever there are com- 
ponents in the system due to undesired wave numbers 
or undesired frequencies. The temporal feedback is im- 
portant both because practical implementations of the 
spatially filtered feedback necessarily involve time delays 
and because spatially filtered feedback alone is sometimes 
not sufficient to stabilize the desired state. 

We emphasize that the control scheme investigated 
here is a plausible candidate for implementation in exper- 
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imental systems. As in Ref. a form of time-delay feed- 
back is used which includes previous comparisons made 
between the state and a time-delayed version in a way 
that is easy to implement because the feedback signal can 
be generated recursively. This is a generalization of 
a low-dimensional control technique known as extended 
time-delay autosychronization, jsj which is in turn a gen- 
eralization of a scheme introduced by Pyragas and has 
been demonstrated to work in low-dimensional electronic 
systems, js] 

Optical systems are of particular interest with respect 
to control both because they offer excellent laboratories 
for testing theoretical ideas and because important tech- 
nological problems associated with them may be solved 
through control techniques. Wide aperture semiconduc- 
tor lasers, with their compact size and very large gain, 
are ideal candidates for high brightness coherent steer- 
able laser sources (spatially and temporally coherent). 
However, the pronounced asymmetry in their gain and 
refractive index spectra leads to a very strong nonlin- 
ear amplitude-phase coupling in the laser field. Con- 
sequently, wide aperture semiconductor lasers display 
uncontrolled dynamic intensity filamentation (random 
beam steering) even immediately beyond lasing thresh- 
old. This behavior persists and becomes even more com- 
plicated at higher current pumping levels. Moreover 
the timescales involved in the laser dynamics lie in the 
nanosecond to picosecond regime, ruling out any algo- 
rithm which requires indirect intervention in order to es- 
tablish control. Time delay feedback with real-time filter- 
ing in space and time is a natural candidate for all-optical 
control of these systems. 

To illustrate the power of spatially filtered, time de- 
lay feedback, we analyze the important example of the 
laser Swift-Hohenberg equations which approximately 
describe the dynamics of the optical field of a wide aper- 
ture semiconductor laser with one transverse dimension. 
Results of linear stability analyses are used to guide the 
choice of control parameters for numerical simulations, 
which reveal that the controlled state may be attained 
even when the initial conditions are far from the lin- 
ear regime. Also, because semiconductor edge-emitting 
lasers typically run on multiple longitudinal modes, we 
study the case of a two-longitudinal mode laser Swift- 
Hohenberg model. We find that the presence of a sec- 
ond mode introduces new features relevant to the control 
scheme, but that control is still possible. 

Our primary motivation for studying this particular 
system is that the feedback signal of interest can be pro- 
duced using a Fabry-Perot interferometer containing a 
spatial filter. The time delay in the feedback scheme 
corresponds to the round trip transit time in the inter- 
ferometer. The spatial filtering can be accomplished by 
a focusing lens, whose focal plane contains the far-field 
fluctuating output of the laser. Placing a suitable aper- 
ture in the focal plane of the lens acts as a narrow band 
spatial filter. One example of a suitable arrangement for 
generating the desired feedback is shown in Fig. |l|. Our 



results therefore suggest a feasible approach to the sup- 
pression of unwanted spatiotemporal fiuctuations in real 
laser systems. 

One of the main results of the present paper is that 
the addition of spatial filtering to time-delay control (us- 
ing An rather than An in Eqn. |^) produces a highly ro- 
bust control scheme. In the context of the laser equa- 
tions discussed below, linear stability analysis of the in- 
finite system shows that in general both the spatial fil- 
tering and the time delay are important components of 
the scheme. For the models we study, numerical results 
also show that when a state is stable with feedback it is 
also highly attracting, so that linear stability analysis is 
predictive even far from the linear regime. Simulations 
of the model equations show that the feedback is able 
to direct the system towards the desired state from a 
distant initial condition, and that spatial filtering is the 
dominant mechanism responsible for this behavior. 

Before proceeding to the detailed analysis of our 
scheme, we mention some related investigations. First, 
results of numerical simulations of the application of our 
scheme with i? = to the control of traveling wave solu- 
tions in the single longitudinal mode semiconductor laser 
Swift-Hohenberg equations (discussed below), have been 
reported elsewhere. Q There, the spatial filtering was 
shown to be extremely robust, rapidly suppressing the 
broadband noisy spatial spectrum of the free running 
unstable laser, in favor of the filtered transverse trav- 
eling wave mode. The subsequent evolution towards a 
controlled state was observed to depend sensitively on 
whether the system is infinitely extended (idealized) or 
pumped over a finite transverse cross-section. In the for- 
mer case, the system evolves to a pure nonlinear trav- 
eling wave mode of the isolated laser system, although 
at higher stress the system typically spent time in a 
metastable dynamical state before reaching the desired 
traveling wave. For finite transverse pumping, the sys- 
tem evolves to a mixed traveling wave solution of the joint 
laser-feedback system and remains there. In this latter 
case, a finite amount of energy remains in the feedback 
loop. 

Second, Lu et al. 0| considered feedback constructed 
by combining comparisons of the current state of the sys- 
tem both to a version delayed by the temporal period of 
the target state and to versions shifted by the spatial pe- 
riod of the target state. Numerical integration showed 
that the control can direct the system towards, and sta- 
bilize, a pattern in a transversely extended laser model, 
but the method does not appear to be a good candidate 
for all-optical implementation. 

Finally, Battogtokh et al. ||] considered the effect of 
feeding back a time-delayed signal constructed from the 
global average of a dynamical field, showing that stable 
uniform oscillatory states of the system with feedback ex- 
ist for some choices of the delay time. In the scheme they 
analyze, the control signal does not represent a small per- 
turbation and the delay time is not tuned to the period 
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of the desired orbit. 
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FIG. 1. (a) Geometry of a wide aperture laser. Ly is as- 
sumed to be sufficiently small that only one mode dominates 
in the y-direction. The large width gives rise to many 
transverse modes. The linear stability analysis in this paper 
is valid for arbitrarily large Lx- The size of determines 
the number of longitudinal modes that are relevant for the 
dynamics, (b) One possible schematic implementation of the 
feedback mechanism studied in this paper. The feedback sig- 
nal is the field reflected from the front of the laser cavity in 
the presence of an additional reflecting interface (labeled "In- 
terferometer Boundary"). The spatial filtering is performed 
by the two lenses with a filter placed at the focal plane of 
each. 



II. SINGLE LONGITUDINAL MODE LASER 
SWIFT-HOHENBERG EQUATIONS 

We now treat the specific example of a recently derived 
model of the transversely extended semiconductor laser, 
the semiconductor laser Swift-Hohenberg equations, 
extending the results of Ref . [01 . The model assumes the 



cavity geometry shown in Fig. Ha) with Ly and both 
small enough that the dynamics is dominated by a single 
mode in the y and z directions, but large. Denoting 
the x-dependent envelope of the electric field by the com- 
plex field ip and the carrier density by the real field n, 
the equations are 



(fj -I- 1)94-0 = cr{r - l)ip + iaV^V - icr^^V' 



dtn 



+ 

— cr(l + iajn^p + 



(2) 



(3) 



where a and b are the decay rates of the electric field and 
population inversion respectively, normalized to the de- 
cay rate of the polarization, r is the scaled pump rate, a 



is a scaled diffusion constant, is the detuning between 
the atomic and carrier frequencies, and a < is a non- 
linear amplitude-phase coupling. (All the coefficients in 
the equation are real.) 

The model is similar to the laser Swift-Hohenberg 
equations for two- level lasers. The key difference 
in the semiconductor equations is the explicit inclusion 
of the a term which derives from the strong assymmetry 
in the semiconductor optical gain and refractive index 
spectra | pT| , p^ . Other terms arising in the semiconduc- 
tor version due to spectral hole burning in the carrier 
distributions do not influence the results discussed here. 

With control turned off (e^ = 0), Eqns. (||) and (||) 
have traveling wave solutions that are always unstable 
iiol: 



ipkix, t) = pexp i(kx — Lot + i 



(4) 
(5) 



where (j) is an arbitrary phase that will henceforth be 
assumed to be zero. 



9^ = b 



r - 1 



l + (J 



LO = 



aO, -\- ak"^ + aap'^/b 



(6) 



(7) 



Note that p is real and that the traveling wave solution 
ceases to exist when the right-hand side of Eqn. (|^) is 
negative. 

Eqn. (||) for the envelope of the electric field contains a 
time-delay and spatial filtering control term of the form 
of Eqn. (^, with the time delay t set to 2tt/u;, the period 
of the desired traveling wave state. The insertion of the 
feedback as simply an additive term in this equation is 
an approximation of the real effect of the feedback, which 
actually consists of an electric field applied at the front 
face of the laser cavity due to reflections from the ele- 
ments shown in Fig. 0(b) . The gain 7 should be thought 
of as a phenomenological parameter that characterizes 
the effect of this boundary term on the longitudinal mode 
in question. The optimal choice of the filter function / 
is not immediately clear. For now we take f{q) to be a 
gaussian of width F; i.e., f{q) — exp[— q^/F^]. As ex- 
plained below, the results are not sensitive to the precise 
choice of F. The case of a square filter function is also 
discussed briefly below. 

To perform the linear stability analysis we write 
ij{x,t) = [1 + B{x,t)]i:k and n{x,t) = [1 + C{x,t)]nk, 
and arrive at the following linearized equations in the 
vicinity of a traveling wave solution: 



(a + l)dtB 



-{2ak + 4iakan ~ Aia^ak^)VB 
+ {ia - 2aan + &da^k'^)V'^B 
-Aiaa^kV^B - aa^V^B 

-(1 + ia)—aC + e_B 


dtC = b{B + B* -C), 



(8) 
(9) 
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where (j = cr/(l + cr)^ 



(12) 
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FIG. 2. (a) Growth rates of perturbations of an uncon- 
trolled (7 — 0) traveling wave solution of Eqns. (8,9) with 
k — 5 and r — 1.5. and parameter values a = .1, n — .001, 
a — .01, 6 = .01, and a = —5. For each q, Re[A] > implies 
exponential growth of the perturbation, (b) Stable region in 
the q, 7 plane of the same traveling wave with control. The 
solid line is the stability boundary for R — 0, t = 27r/tj, and 
r — .25. Note that all modes are stable for 7 < —0.1. The 
dashed line (see inset) corresponds to the (unphysical) case 
of r = 0, i? = 0. 

Following Ref. we obtain a linear system of ordi- 
nary differential equations for the Fourier modes of the 
perturbation. Letting ^ = {Bq, B*, Cg), the 3D vector of 
Fourier amplitudes at wave number q, the equations can 
be written in the general form 



(10) 



where e is given by the expression in Eqn. (n), J is ob- 
tained from the coefficients of Eqns. (||) and (9|), and M 
is determined by which variables form the control signal 
and how the control signal enters the equations. In the 
present case, 



M = 




(11) 



The factor by which a given eigenmode of the per- 
turbation grows during one period of the evolution of the 
controlled system is called a Floquet multiplier. The time 
delay in the control term requires that the initial condi- 
tions for the evolution must specify the behavior over an 
entire continuous time interval of one period, so each spa- 
tial Fourier mode has an infinite number of eigenmodes 
and Floquet multipliers, /ij. Letting ^'^■'^ represent the 
j*^ eigenmode, we have, by definition. 



The set of Floquet multipliers for a perturbation with 
wave number q determine that perturbation's linear sta- 
bility; if one or more multiplier has \fj,j\ > 1, the pertur- 
bation is unstable. 
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FIG. 3. Evolution of the single mode system with 
a, b, (J, a, r, k as in Fig. 2a. , and 7 = .2. (a) A spacetime 
plot of the phase of the field ^l} in the uncontrolled system 
showing chaotic fluctuations, (b) The growth of the magni- 
tude of the desired mode as a function of time, (c) A space- 
time plot of the phase of ijj in the controlled system for the 
same run as shown in (b). The lower region shows the dy- 
namics when the system is first turned on. After a short 
period during which fluctuations grow rapidly, the feedback 
suppresses all wave numbers other than the desired one. After 
a transient time of approximately 500 periods of the desired 
orbit, the system settles into the traveling wave state. The up- 
per portion shows the system as as it approaches the desired 
state, which would appear as straight bands on this plot. 

Dropping the subscript j and evaluating the geometric 
sum in e, Eqn. (Hn) may be written 



(13) 



The values of the Floquet multipliers are determined 
by requiring consistency between this equation and the 
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defining relation of Eqn. (|12|). We obtain the following 
characteristic equation for this modified eigenvalue prob- 
lem: 



J+7 



l--fi/2(q)n- 



M 



0, (14) 



where the exponential represents the operator that ad- 
vances the linear system by one period r. As discussed 
in Refs. jl^,^ , one can perform a numerical winding num- 
ber calculation oi g{^~^) around the unit circle to obtain 
the number of roots satisfying < 1. Since there are 

no poles in the unit disk, the system is linearly stable if 
and only if this winding number vanishes. 

Results from the linear stability analysis predict that 
our control technique successfully stabilizes all traveling 
wave solutions in the single longitudinal mode model. 
We present detailed results for a single traveling wave 
solution, for k — 5 and r = 1.5, which is typical of all 
traveling waves we have studied. (Values of the other 
parameters are given in the caption.) 

Fig. Ija shows the growth rates of the modes of the 
uncontrolled system, which are obtained by finding the 
eigenvalues of J in Eqn. |l^. There is one unstable mode 
for perturbation wave numbers between zero and ^ 12. 
With control, using i? = 0, we find that the traveling 
wave state is stable for 7 sufficiently negative. The solid 
line in Fig. |^ indicates the boundary between which per- 
turbation wave numbers are stable or unstable at a given 
7. The controlled traveling wave is stable at values of 
7 for which all wave numbers are stable, i.e., where the 
shaded region contains an entire horizontal line. For all 
traveling waves in this model, there is a minimum I7I for 
which the state is stable. In the case shown in Fig. ^3, 
this occurs at 7 ~ —.1. In this model, there is no lower 
boundary to the stable region for traveling waves. 

We find for this system that spatial filtering alone 
would be sufficient to stabilize the traveling wave. The 
stability boundary obtained with t — and i? = 0, the 
dashed line in Fig. ||b, is nearly identical to the one with 
T = In ju) at large q, but the time delay clearly has a 
significant effect at q near zero. It is also important to 
note that implementation of a spatial filter with no time 
delay is not possible in fast optical systems. The result 
that the introduction of a time delay of one period does 
not destroy the stability in the case of the gaussian filter 
is therefore significant. 

In general, a given wave number perturbation can be 
stabilized either by the time delay feedback with no spa- 
tial filtering or by the spatial filtering with no time delay. 
In each case, however, there may be small bands of wave 
numbers for which one or the other method fails. A given 
spatial filter fails near g = if there exist perturbations 
which are sufficiently unstable (or if /(g) is sufficiently 
close to unity). For g's at which /(g) = 1, as occurs 
for a finite band in the step function case, the spatial 
filtering has no effect on the stability. The time delay 
feedback alone fails for wave numbers whose frequency 



of oscillation is sufficiently close to an integer multiple of 
the frequency of the desired traveling wave. Combining 
the spatial filter and the time delay renders the system 
stable at all q. 

The time delay is a crucial component for stability in 
the two mode system discussed below. In the single mode 
system, it may also play an important role if /(g) is cho- 
sen to be a step function rather than a gaussian. 

The predictions of several stability diagrams similar to 
Fig. 1^ have been checked in detail by numerical simula- 
tion. The numerics show that the traveling wave states 
are stabilized with values of 7 predicted by the linear 
analysis, and that instabilities occur at the wave number 
predicted when I7I is too small. 

An important question is whether the linear stability 
analysis is predictive of the behavior of the system even 
for initial conditions that are not in the linear regime. 
Numerical integration of the model equations show that 
the spatially filtered feedback is particularly effective in 
directing the system to the desired state. As illustrated 
in Fig. 0, for parameters corresponding to a linearly sta- 
ble controlled state, the system is attracted to the de- 
sired state from a typical initial condition. Though it 
is difficult to display the full behavior during the long 
transient, an investigation of the details reveals that, be- 
ginning from low amplitude noise of the type that would 
be expected when the laser is first turned on, the system, 
depending on the parameter regime, may pass through 
several nearly stable states with the desired wave num- 
ber, but the incorrect frequency, before finally settling 
on the one with the desired frequency. Preliminary in- 
vestigations of systems with time delay feedback alone 
indicate that more complicated behavior occurs beyond 
the linear regime. 



III. TWO LONGITUDINAL MODE LASER 
SWIFT-HOHENBERG EQUATIONS 

Semiconductor lasers of practical interest generally op- 
erate in regimes where many longitudinal modes may be 
active. To begin to understand the possible effects of 
multiple longitudinal modes, we study a two mode model. 
This model is a straightforward generalization of the two 
level, one mode model derived in Ref. |l^ to the situation 
in which two longitudinal modes, with mode separation 
2A, dominate the dynamics. |14| With the addition of 
the semiconductor a term discussed above, this model 
reads, 



(cr -I- V)dti)\ = o'('' ^ l)V-'i + iaV^Vi + iA^/>i — iafl-ipi 



{1 + ar 

-(7(1 + ia)nipi - a{l + ia)r]ip2 + e^, (15) 
(cr + l)(?tV'2 = o-(r - 1)V'2 + «aV^V'2 - «AV'2 - io-ilV'2 



{1 + ay 



5 



dtn = 



-bn + \ipi^ 



ia)'mjj2 — 



|2 



ia)ri*%l}i + (16) 
(17) 
(18) 



Note that the same control term = e^^ + appears 
in the both if^i and "02 equations with equal magnitude. 
This simple way to model the effect of the reinjection of 
the reflected field into the laser cavity is used here for 
convenience. The present model is intended only to dis- 
play the new qualitative features that arise when more 
than one mode is relevant. 
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FIG. 4. Dispersion curves for (unstable) traveling wave so- 
lutions in the two mode model. The solid (dashed) line rep- 
resents solutions in which only the favored (unfavored) mode 
is excited. The parameters are 5 = .05, a = .1, n = .001, 
a — .01, b — .01, and a — —5. 

We are interested in the solution in which one longi- 
tudinal mode supports a traveling wave and the other is 
inactive: 



n 

i>2 



ipk 

nk 
0, 



; pexpi(fca; — wt), 



77 = 0, 



where 



crfi -I- aA:^ - A + aap^/b 



1 



(19) 
(20) 
(21) 
(22) 



(23) 



(24) 



The complementary solution is obtained by interchang- 
ing the subscripts of the fields and taking A — > —A in 
the expresions for p and uj. Taking ■01 = (1 + i?)'0fe, 
1^2 = Dipk, n = (1 + C)nk, and rj = E, we obtain the 
following linear equations for the small fields B, D, C, E, 



(a + l)dtB = -{2ak + Aiakdn - 'iia^dk^)VB 

+ {ia-2aan + 6aa'^k'^)W^B (25) 
(cr + l)dtD = [zaV^ - 2ik A + aV^ f 



-iak^ 



CT(r2 + A - ak^ f] D 



dtC 
dtE 



-<j{l + ia)E* 
b{B + B* -C) 
-2bE + p^D* 



(26) 
(27) 
(28) 
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FIG. 5. (a) Real parts of the eigenvalues of an uncontrolled 
(7 = 0) solution in which 1P2 is a traveling wave with fc = 8 
and r = 1.5, and tpi is zero everywhere. Parameters are the 
same as in Fig. 4. (b) Stable region of the same solution with 
control with parameters R = .5, t = 2n/u>, and F — .25. The 
traveling wave is stable at all q for —0.13 < 7 < —0.05. 

Fourier transforming, we again obtain a general ex- 
pression for the behavior of small differences of a pertur- 
bation wave number from the controlled state. Letting 
Qx be the wave number in the transverse direction and 
e = {B, B*,D, D*, C, E, E*), we have 



M- 



(29) 



where here J is the matrix of coefficients obtained di- 
rectly from Eqns. ( p5| - |2^ ) and 



/ 1 1 \ 

10 10 

1 1 



M 



10 10 



Voooooooy 



(30) 



As in the case of the single longitudinal mode laser, a 
condition of the form of Eqn. 14 defines the Floquet mul- 
tiphers of the system. 



6 



a 



0.00 



-0.05 



-0.10 



unstable 




0.00 0.25 

Perturbation Wavenumber q 



0.50 



FIG. 6. A detail from the left edge of the upper bound- 
ary of the domain of control in Fig. 5(b) shows the effect of 
varying R. In the case of R = 0, as with no time delay, the 
traveling wave solution cannot be controlled due to instabil- 
ities at small q. Larger values of R do yield stable solutions 
for sufficiently negative 7. 

We now describe the results of the linear stability anal- 
ysis of the two mode model. As in the single mode 
model, each mode is always unstable to transverse fluc- 
tuations, but in the two mode model it is possible for one 
mode to be unstable to the growth of the other as well. 
A straightforward stability analysis of the uncontrolled 
equations shows that for all parameter choices both ipi 
or -02 are marginally stable against transverse fluctua- 
tions at qx = 0, but only one of the modes is always 
stable against growth of the other mode. Which mode is 
which depends upon the choice of the mode separation, 
the desired wave number, and other parameters in the 
model. We will refer to a mode that is stable (unstable) 
against growth of other longitudinal modes at = as 
"favored" ( "unfavored" ) . 

The dispersion curves for transverse waves in the two 
modes have nearly the same functional form, but are dis- 
placed relative to each other approximately by the mode 
spacing, 2A. (See, for example. Fig. ^) By choosing the 
wave number for the spatial filter, one selects one trav- 
eling wave state from each of the two dispersion curves. 
Because the frequencies of these two states are different, 
one can choose the time delay so as to suppress fluctu- 
ations at the frequency of the undesired mode. Thus it 
is plausible to suppose that the combination of the time 
delay and the spatial filter capable of stabilizing either 
of the two longitudinal modes. We will focus on the sta- 
bilization of an unfavored mode, both because it would 
appear to be the more difficult case and because it may 
be a better representation of the situation that arises in 
multimode systems. 

Fig. H illustrates the stabilization of the unfavored 
mode (ipi) with the parameters listed in the caption. 
Fig. ^a, shows the stability curves for the uncontrolled 
system, clearly indicating the instability at g = that 
makes this case qualitatively different from the single 



mode case discussed above. The spatial filter compo- 
nent in our control scheme is insensitive to instabilities 
at or near q — because the filter must pass components 
of both ipi and ^2 with this wave number. The width 
of the function f{q) used for the spatial filter will de- 
termine the range of q which are passed. As a result of 
the ineffectiveness of the spatial filter over this range of 
perturbation wave number, the temporal component of 
the control scheme must be relied on to stabilize these 
perturbations. 



o 

T3 
CD 

O 
> 
TO 

H — 

ZD 




Favored 



ro 
o 

T3 

CD 

o 
> 

03 



0^ 

200 400 600 800 
Time (Periods of Traveling Wave) 
(a) 




-6 -4 -2 2 4 6 
X (Transverse a.u.) 

(b) 

FIG. 7. Evolution of two mode system with a, b, Q, r, a, a, 
as in Fig. 4. with A = 0.1, r = 1.5, k = S, R ^ 0.5, and 
7 = .1. (a) A spacetime plot of the phase of in the con- 
trolled system. The lower region shows the dynamics when 
the system is first turned on. After a transient time of ap- 
proximately 200 periods of the desired orbit, the system set- 
tles into the traveling wave state, (b) The magnitude of the 
favored mode (dashed line) and unfavored mode (solid line) 
as a function of time for the same run as shown in (b). Note 
the expanded scale at the right for the favored mode. 

The stability diagrams of Fig. ||b and Fig. ^ demon- 
strates that the time-delay control is effective in control- 
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ling the range of perturbation wave numbers that are not 
stabihzed by the spatial filter. Fig. ^ shows that with 
both time-delay control (here with R = 0.5) and the 
spatial filter (with T — 0.25) there is a range of 7 that 
stabilizes the traveling wave solution. Fig. M shows that 
when time-delay control is not present, and also when 
R is too small, there is no range of 7 that stabilizes the 
traveling wave solution at all wavenumbers. 

A new feature that appears in the two mode model, 
and is shown in Fig. is the lower boundary of the sta- 
ble domain, whose origin lies in the off-diagonal elements 
of M. When the system is not exactly on the desired 
orbit, there is a finite amount of feedback generated. Be- 
cause the desired mode has a much larger average mag- 
nitude than the other mode, the feedback signal is domi- 
nated by effects from the desired mode. This feedback is 
necessary to control the desired mode, but it also affects 
the other mode. When the magnitude of this feedback 
becomes too large, as it must when I7I is increased, these 
unwanted perturbations to the undesired mode cause the 
state to go unstable. 

The position of the lower boundary of the domain of 
control (Fig. ^) is important because it determines the 
range of gain that can be used to obtain control. If that 
range is very small, it may be difficult to find an ap- 
propriate 7 in an experiment. Even worse, if the lower 
boundary becomes so high that part of it reaches the 
lowest point of the upper boundary, there is no 7 which 
can control the system. We find that the position of the 
lower boundary is affected by several parameters. The 
lower boundary is raised when the pump rate r is raised 
and when the wave number k is lowered. The mode sep- 
aration, A, also plays an important role in the location 
of the lower boundary. For larger A, the lower bound- 
ary is pushed down. In a system in which 7 is the only 
adjustable parameter (r. A, and R fixed), we find that 
traveling waves with wavenumbers in a finite continuous 
band can be stabilized. The high-fc boundary of the band 
is determined by the condition that traveling waves ex- 
ist (that p must be real), and the low-fc boundary is the 
point at which there ceases to be a 7 that can control 
perturbations at all wave numbers. 

As in the single mode model, numerical simulations 
confirm the predictions of the linear analysis and show 
that the traveling wave state can be obtained starting 
from a distant initial condition. Fig. |^ shows the emer- 
gence of the desired traveling wave from a low amplitude, 
noisy initial condition. After an initial transient, the sys- 
tem clearly settles into the desired pure traveling wave. 

We have also observed the behavior of the system when 
R is chosen too small. Although the only unstable modes 
in this case are very close to q = 0, we find that their 
growth completely destroys the traveling wave. The sys- 
tem does not merely develop long wavelength modula- 
tions of the desired wave. We therefore conclude that 
both the temporal and spatial aspects of the feedback 
signal we have analyzed play essential roles in the suc- 
cess of the scheme. 



IV. CONCLUSIONS 

Our study of the dynamics of laser Swift- Hohcnberg 
equations with time-delayed, spatially filtered feedback 
strongly suggests that stable lasing at a single transverse 
wave number in wide aperture lasers is possible. A fu- 
ture publication will report on studies of a more realistic 
model of field and carrier dynamics in a semiconductor 
system of the type shown in Fig. 0, where preliminary 
results are encouraging. Though there are several non- 
trivial experimental issues associated with the fabrication 
of such a device, we believe that this is a promising di- 
rection for research and development. 

We have presented a theoretical approach to the anal- 
ysis of this sort of feedback that appears to capture the 
relevant features of the dynamics. The linear stability 
analysis presented here is a straightforward extension of 
some of previous work on stabilizing traveling waves in 
the complex Ginzburg-Landau equation. In the present 
case, however, the desired state seems to be a global at- 
tractor, which gives us considerably more confidence in 
its potential for practical implementation. 

Finally, we would like to emphasize that the general 
method of applying time-delay feedback combined with 
spatial filtering is a powerful technique that might be 
adapted to many other types of physical systems. Its pri- 
mary advantage is that the desired traveling wave state 
need not be available in some external form for construc- 
tion of the feedback signal. It is particularly suitable for 
optical systems, however, where the necessary manipula- 
tions of the signal can be performed with standard optical 
elements. 
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